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A search for long-lived neutral particles (A*'") which decay into at least one muon has been 
performed using an instrumented decay channel at the E815 (NuTeV) experiment at Fermilab. The 
decay channel was composed of helium bags interspersed with drift chambers, and was used in 
conjunction with the NuTeV neutrino detector to search for A*''^ decays. The data were examined 
for particles decaying into the muonic final states fijj,, fj.e, and ij.n. Three fj.^ events were observed 
over an expected background of 0.040 ± 0.009 events; no events were observed in the other modes. 
Although the observed events share some characteristics with neutrino interactions, the observed 
rate is a factor of 75 greater than expected. No Standard Model process appears to be consistent 
with this observation. 



I. INTRODUCTION 

In various extensions to the Standard Model, new parti- 
cles exist which have reduced couplings to normal quarks 
and leptons. These new particles may have zero electric 
charge, long lifetimes, and small interaction rates with 
normal matter. We shall refer to these as iV° particles 
in the following text. Examples of such N'^ particles 
include neutral heavy leptons (NHLs) or heavy sterile 
neutrinos [P-pl and neutral supersymmetric particles B 
such as neutralinos and sneutrinos. The N*^ particles can 
be produced either by pair production in hadronic inter- 
actions or via weak decays of mesons through mixing 
with standard neutrinos. The decays of the N'^ to nor- 
mal hadrons and/or leptons can proceed through weak 
decays with mixing, or via i?-parity violating supersym- 
metric processes. 

High energy neutrino beamlines are ideal places to pro- 
duce N'^ particles, since very large numbers of protons 
interact in these beamlines. A^'^'s may be produced via 
a number of mechanisms, including primary interactions 
of the protons either in the target or the beam dump, 
through prompt decays of charmed or bottom mesons, 
by decays of pions or kaons in the decay region, or in 



neutrino interactions in the shielding downstream of the 
decay region. A particle detector placed downstream of 
this sort of beamline (i.e., in the neutrino beam itself) 
can be used to search for N° decays. 

We report here the results of a search using Fermilab's 
E815 (NuTeV) detector for N*^ particles in the mass re- 
gion above 2.2 GeV/c^ which decay into final states with 
at least one muon and one other charged particle. For 
the search described here, the NuTeV neutrino beamline 
was used in conjunction with a low mass decay detector 
called the decay channel. 

NuTeV has previously reported results of searches for 
A^^'s in the mass region between 0.3 to 3.0 GeV/c^ with 
at least one final state muon H , and in the mass region 
below 0.3 GeV/c^ for decays to electrons Q. The 0.3 
to 3.0 GeV/c^ study addressed NHLs that could be pro- 
duced in the decay of K and D mesons, whose hadronic 
production rate is known M. This mass region also has 
low background from deep inelastic neutrino events in 
the decay channel. The low mass (< 0.3 GeV/c^) study 
was pursued mainly to address the KARMEN timing 
anomaly B , which has been interpreted as a N^ particle 
with a mass equal to 33.9 MeV/c^. 

The search for events with masses above 2.2 GeV/c^ 
(which we shall refer to as "the high mass region" ) is dif- 
ferent from the previous searches in two respects. First, 



the backgrounds from neutrino interactions are much 
higher in the high mass range. Reducing this poten- 
tial background source required tighter selection criteria 
than the ones used in the previous analyses. Second, the 
production mechanisms are quite different from those for 
lower mass NHLs. NHLs with masses above 1.8 GeV /c? 
arise either from decays of B mesons or from neutrino 
interactions in the neutrino beam shielding ( "berm" ) up- 
stream of the detector |^] . The B meson production cross- 
section at 800 GeV/c is much smaller and less well known 
than that for the lighter mesons, and this partially moti- 
vated separating the high mass search from the previous 
ones. Recently, though, new data from Fermilab E771 
have indicated a larger B production cross-section, mak- 
ing a search for NHLs from this source feasible ||lO(| . 



II. THE BEAMLINE AND DETECTOR 

During the 1997 fixed-target run at Fermilab, NuTeV 
received 2.54 x 10^^ 800 GeV/c protons with the detector 
configured for this search. The proton beam was incident 
on a one-interaction-length beryllium oxide target at a 
targeting angle of 7.8 mr with respect to the detector. 
A sign-selected quadrupole train (SSQT) jlll focused ei- 
ther positive (for 1.13 x 10^^ protons) or negative (for 
1.41 X 10^* protons) secondary tt and K mesons into a 440 
m evacuated decay region pointed towards the NuTeV de- 
cay channel and neutrino detector hall. Surviving neu- 
trinos (and possibly also A^^'s) traversed ~850 meters of 
earth-berm shielding before reaching the NuTeV decay 
channel. 

The decay channel region (Figure [|), located 1.4 km 
downstream of the production target, was designed to 
contain minimal material (in order to suppress neutrino 
interactions) and to have tracking sufficient to isolate 
two-track decays of neutral particles. A 4.6 m x 4.6 m 
double array of plastic scintillation counters vetoed 
charged particles entering from upstream of the decay 
channel. If two counters in a back-front coincidence fired, 
the event was vetoed. Timing resolution for the veto sys- 
tem was 3.8 ns. A NuTeV testbeam chamber was po- 
sitioned immediately downstream of the veto wall and 
offset from the center of the decay channel; this chamber 
was not used in this analysis. The channel itself mea- 
sured 34 m in length and was interspersed with 3 m x 3 m 
argon-ethane drift chambers positioned at 14.5 m, 24 m, 
and 34 m downstream of the veto array in stations of 1, 1, 
and 4 chambers, respectively. From upstream to down- 
stream, the chambers were labeled DK5 through DKl, 
followed by TG43 at the front face of the neutrino detec- 
tor. Chambers DK5 through DK2 and TG43 had a single 
sense wire per cell ( "single- wire" chambers) , resulting in 
a left-right ambiguity in the hit position. This ambiguity 
was reduced in tracking by making use of the fact that 



the positions of the single-wire chambers were staggered 
in X and y. Chamber DKl was of a "three-wire" (two 
sense wires and one field-shaping wire per cell) design, 
which helped resolve any remaining tracking ambiguities. 
DKl was also rotated by 47.7 mr about the beam axis, 
allowing tracks in the x view to be matched with those 
in the y view. The regions between the drift chamber 
stations were occupied by hclium-fiUcd cylindrical plas- 
tic bags 4.6 m in diameter. 

In the offline analysis, tracks were reconstructed from 
drift chamber hits and grouped together to form ver- 
tices. The tracking algorithm took into account multiple 
Coulomb scattering, using a full error matrix for the fit. 
Sets of tracks were grouped as candidates for a vertex 
if their distance of closest approach was less than 12.7 
cm. The vertex position was then determined using a 
constrained fit. Typically, a vertex from a N^ of mass 5 
GeV/c^ would be reconstructed with a resolution of 0.13 
cm in the transverse direction and 7.4 cm longitudinally. 

The Lab E neutrino detector [|l^,|3|, located imme- 
diately downstream of the decay channel, provided fi- 
nal state particle energy measurement and identification. 
This detector consisted of a 690 ton iron-scintillator tar- 
get calorimeter followed by a toroidal niuon spectrome- 
ter. Three-wire argon-ethane drift chambers were posi- 
tioned every 20 cm along the length of the calorimeter, 
and 84 2.5 cm-thick liquid scintillator counters were inter- 
leaved with the steel plates at 10 cm intervals throughout 
the its length. The spectrometer had a 15 kC toroidal 
magnetic field with drift chambers interspersed through- 
out the toroid magnets to provide tracking for the muons. 

Sets of hits in the calorimeter drift chambers were 
linked to tracks found within the decay channel to de- 
termine the particle identification for each track. By an- 
alyzing the distribution of hits in the calorimeter as either 
single, long tracks (consistent with a muon), as a com- 
pact cluster of hits (consistent with an electron shower), 
or as an elongated cluster of hits (consistent with a pion 
shower), the particle type was determined. All charged 
hadrons were reconstructed as pions; "/^tt" samples re- 
ferred to below are understood to include any fip or fiK 
final states as well. The NuTeV calibration beam pro- 
vided electrons, pions, and muons to the Lab E detector; 
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FIG. 1. Schematic of the NuTeV decay channel with ex- 
ample N^'' -^ /iTT decay. 



TABLE I. Misreconstruction rates for Monte Carlo N° de- 
cays 



Generated 
Decay Mode 



Reconstruction mode probability 



MM 



fie 



flTV 






100+°:;^% 



0.0 /» 



o.o: 



% 



f0.3o 
-0.0' 



0.0 

1±1% 



89 ± 2% 
25 ± 3% 



11 ± 2% 

74 ± 3% 



these were used to tune the particle identification algo- 
rithm. Misidentification rates for 7V° decays were deter- 
mined using the GEANT-based Monte Carlo (MC) fll] 
described in Section 0; these rates are listed in Table]|! 

In the case of electrons and pions, pulse height infor- 
mation from the counters was used to determine par- 
ticle energy deposition. The hadronic energy resolu- 
tion of the calorimeter was a/E = (0.024 ± 0.001) © 
(0.874 ± 0.003)/-\/E; the electromagnetic energy resolu- 
tion, a/E = (0.04 ± O.OOl) ® (0.52 ± 0.01)/\/^ (l|. If 
two or more pion- or electron-clusters were present, the 
energy determined from the pulse heights was divided 
according to the number of drift chamber hits in each 
cluster. 

Muon energy determination depended on the topology 
of the track. If the muon track extended into the toroid, 
the spectrometer measurement was used (resolution of 
11%). If the muon stopped in the calorimeter steel, the 
momentum was determined by range (resolution of 310 
MeV). The energy of muons which exited through the 
side of the calorimeter was determined from the track's 
multiple scattering in the steel. For a 50 GeV/c muon, 
the resolution for this method is 42%. 



III. PHILOSOPHY OF THIS ANALYSIS 



were used. In those cases where new cuts were intro- 
duced, demonstration of a strong MC-based motivation 
was required. New members who had not seen the event 
joined the analysis group. Finally, an important aspect 
of the analysis included setting up orthogonal analysis 
regions and comparing Monte Carlo prediction to the re- 
sult. Each of the two previously- published analyses was 
motivated by its own physics goals, but they also repre- 
sent tests in regions complementary to this analysis. 

As part of the analysis philosophy, once the analysis re- 
gion was selected based on the Monte Carlo criteria, the 
collaboration agreed to show any events which were ob- 
served. However, the interpretation of the events might 
change after the analysis region was examined, upon fur- 
ther investigation. 



IV. EVENT SELECTION 

Event selection criteria were developed to minimize 
known backgrounds while maintaining efficiency for a 
possible 7V° signal. 

As an example of what might be observed in the de- 
cay channel. Figure g shows a Monte Carlo simulated 
event of a 5 GeV/c^ N^ -^ jjl^v on the NuTeV event dis- 
play. The beam enters the decay channel from the left of 
the figure. The decay channel chambers (DK5 through 
TG43) appear sequentially from left to right. The ver- 
tex of the simulated decay is immediately downstream of 
DK4. Hits in chambers are indicated by crosses; two hits 
per track are shown in each chamber because of the left- 
right ambiguity in the single-wire chambers. The muons 
can be seen to penetrate the calorimeter steel and then 
bend in the toroidal magnetic field at the far right. Pulse 
height information from the counters is shown by the his- 
togram above the calorimeter region. 



In analyses such as this, there is a real concern that 
events may be eliminated or isolated through an uninten- 
tional bias of the people involved in the analysis. The so- 
lution adopted by many collaborations is that of a "closed 
box" analysis, in which there is no direct access to the 
signal region until the end of the analysis. This procedure 
was a philosophical goal of this search. However, before 
this analysis, during the early development of the recon- 
struction software for the decay channel, one candidate 
decay channel event with two muon tracks was observed. 
This event was studied in detail and ascertained to have 
a mass greater than 2.2 GeV/c^. 

Because of the observation of this event, the NuTeV 
collaboration went to considerable effort to minimize 
bias. Investigations of data events with high mass were 
stopped until the MC background studies (described be- 
low) were completed to establish the cuts and require- 
ments. In most cases, cuts set prior to the observation 




FIG. 2. Monte Carlo simulated event for a 5 GeV/c^ 



TABLE II. Overview of reconstruction cuts used in the analysis. 



Event conditions 


No veto, physics trigger, protons in the spill, run period 


Cut cosmic ray tracks 


slope < 0.1 radians 


3-D track reconstruction 


X and y views linked for each track 


Track & vertex reconstruction 


Track pseudo-xVdof< 10 




Vertex pseudo-x^/dof< 10 


Require two-track vertex 


Only two tracks connected to upstream vertex 




only one particle downstream. 


Good particle ID 


Track is within calorimeter acceptance 




Cluster associated with track 




/i Energy > 2.2 GeV, e and tt energies > 10 GeV 


Cut Kl punch-through 


T.Etracks > 12 GeV 


Isolate fiducial volume 


vertex within \x\ < 127, \y\ < 127 cm 




vertex 101.6 cm away from chambers 




or 3 times vertex error (whichever is larger). 


Isolate mass region 


Mt > 2.2 GeV 



In decays with a neutral particle in the final state, 
it is not possible to reconstruct the invariant mass of 
the A^'^. Instead, one must use the "transverse mass," 
rriT = I-PtI + \/Pt + "^yi where Pt is the component 
of the total reconstructed momentum perpendicular to 
the beam direction, and my is the invariant mass of the 
visible particles. In the case of perfect resolution, the 
transverse mass of the event is always lower than the N^ 
mass. When there is no final state neutrino, and hence 
no missing Pt, this expression reduces to the invariant 
mass. 

The main sources of conventional events in the decay 
channel are: 1) deep-inelastic scattering (DIS) of neutri- 
nos or anti-neutrinos in the drift chamber material; 2) 
DIS in the helium; and 3) DIS in the material surround- 
ing the decay channel with a misreconstructed vertex in 
the fiducial volume. The DIS event rate rises with Mt, 
and represents the dominant background for the high 
mass analysis. 

Other sources of background are small compared to 
DIS. Neutral kaons produced in surrounding material 
may enter the channel and decay in the fiducial volume. 
These will appear to have high transverse mass if the 
kaon enters with a large angle. Diffractive tt, K, p, and 
charm production from neutrino interactions may occur 
in either the chambers or the helium. Another interaction 
which can occur in either chambers or helium is low mul- 
tiplicity neutrino-induced resonance production, charac- 
terized by a high-energy forward muon accompanied by a 
low-energy pion track. A number of other possible back- 
ground sources have been found to be negligible because 
they very rarely produce reconstructed vertices in the 



decay channel. These include cosmic ray showers, con- 
versions of photons produced in surrounding material, 
interactions from muons scattered from surrounding ma- 
terial, "leakage" of charged particles from the adjacent 
testbeam line (which occurred only during specific data 
collection periods), an out-of-time neutrino interaction 
overlaid on an in-time interaction, and two coincident 
in-time interactions in the decay channel. These back- 
ground sources were constrained by data as well as in- 
vestigated through Monte Carlo. 

The goal of this stage of the analysis was to create a 
set of cuts which reduced the number of events expected 
from conventional sources to well below one event. Most 
of these cuts were originally developed for the 0.3 to 3.0 
GeV/c^ search. The cuts used to isolate N^ decays fell 
into two broad categories: reconstruction and "clean" 
cuts. 

Reconstruction cuts isolated events with exactly two 
tracks forming a vertex within the decay channel fidu- 
cial volume and having no charged particle identified in 
the upstream veto system. A summary of these cuts ap- 
pears in Table |l[ Both tracks were required to be well- 
reconstructed. A "pseudo-X(rfc/dof" (degrees of free- 
dom) was used|^ to measure track reconstruction quality. 
Tracks were required to have an associated calorimeter 
cluster, with at least one of the tracks identified as a 



^ "Pseudo" means that Gaussian errors on hit positions were 
assumed. 



muon. Cuts were also applied on vertex quality (pseudo- 
X^2./dof < 10, which corresponds to a 95% acceptance 
probability) and transverse position within the detector 
fiducial volume (|x| < 127 cm, \y\ < 127 cm). In or- 
der to remove events which might be due to interactions 
in the chambers, events were cut where the longitudinal 
(z) distance from the vertex position to any drift cham- 
ber was less than either ±101.6 cm or three times the 
vertex error. Events were allowed to have a second ver- 
tex downstream of the first to allow for the possibility 
of events with (5-rays. Cosmic ray tracks, which gener- 
ally have large angles with respect to the beam direction, 
were removed by requiring the slope of each track to be 
less than 100 mr. To ensure accurate particle identifi- 
cation and energy measurement, muons, hadrons, and 
electrons were required to have an energy greater than 
2.2 GeV, 10 GeV, and 10 GeV, respectively This lat- 
ter cut also eliminated low-energy pions associated with 
neutrino-induced resonance production. An additional 
total energy cut of 12 GeV was applied to nfi events to 
remove background from K^ -^ t^jJiv decays with a sub- 
sequent decay of the pion. In order to isolate high mass 
events, a cut of nriT > 2.2 GeV/c^ was applied. 

"Clean event" cuts were applied to reduce the deep- 
inelastic neutrino scattering backgrounds. These cuts, 
specific to this analysis, were motivated by the observa- 
tion that DIS events typically have large track multiplici- 
ties with many drift chamber hits and extra, unassociated 
clusters in the calorimeter. A summary of the character- 
istics associated with DIS events wh ich pass the recon- 
struction cuts is given in Table III. The "clean" cuts 
included the following requirements: 1) three or fewer 
tracks in any one view, 2) three or fewer drift chamber 
hits in any view of the first chamber downstream of the 
vertex, 3) at least one view (x or y) with fewer than 
eight drift chamber hits total in the first two chambers 
downstream of the vertex, 4) no energy clusters in the 
calorimeter not associated with tracks, and 5) no tracks 
identified as electrons with missing hits in either view of 
the first two chambers downstream of the vertex. Re- 
quirements 1-3 remove events with high multiplicities; 
requirement 4 removes events where a neutral particle 
deposits energy in the calorimeter; and requirement 5 



was used to reduce events with photon conversions in 
the downstream chambers which could be misidentified 
as electrons. 



V. BACKGROUND ESTIMATION USING 
MONTE CARLO 

Detailed Monte Carlo simulations of both physics pro- 
cesses and detector effects were used to quantify the 
background from neutrino interactions after cuts. In- 
put to the simulation was provided from several event 
generators. The LEPTO/Jetset Monte Carlo program 
was used to simulate DIS events jisll . This simulation 
used CCFR parton distributions ||16[1 , included the cor- 
rect A-dependence [O], and generated DIS events from 
Q^ > 0.1 GeV^ and W > 2 GeV. Resonance and con- 
tinuum production, simulated using the calculations of 
Belusevic and Rein Q, allowed us to extend the Monte 
Carlo into the \ow-W region. Diffractive production 
was calculated using Vector Meson Dominance (VMD) 
and Partially-Conserved Axial Current (PCAC) models 
normalized to a previous measurement with the NuTeV 
calorimeter data [E9l. 

The event generators fed a GEANT-based @ detec- 
tor simulation that produced hit-level simulations of raw 
data. Cell-by-cell inefficiencies and dead regions due to 
internal chamber supports were included. To simulate 
noise and accidental activity in the detector, decay chan- 
nel hits taken from in-time downstream calorimeter neu- 
trino events were overlaid on the GEANT events. Monte 
Carlo events were processed using the same analysis rou- 
tines used for the data. 

Background calculations were normalized to the data 
using charged-current DIS interactions in the chambers. 
Events in this sample were required to pass the following 
five "normalization cuts" : a vertex within the transverse 
fiducial volume (|a;| < 127 cm, |?;| < 127 cm); a z ver- 
tex within 76.2 cm of a drift chamber; no coincidences 
within ±50 ns of the trigger in the upstream veto system; 
>1 GeV energy deposit in the front of the calorimeter; 
and one toroid-analyzed muon matched to a decay chan- 
nel track. The Monte Carlo was normalized to match 



TABLE III. Identification of events with exiting tracks, neutral particles, and photon conversions; these char- 
acterize DIS events and are removed by the "clean" cuts. 



Extra tracks 


>3 tracks in either view 


Exiting tracks 


>1 extra hit in either view in the chamber just downstream of vertex 
or >7 hits in each view in the first 2 chambers downstream of the vertex. 


Neutral particles 


> 1 cluster (s) in calorimeter without associated track. 


Photon conversion 


Electron PID with no hits in chamber immediately downstream of vertex. 
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FIG. 3. Distribution of number of reconstructed decay 
channel particles comparing data (points) and Monte Carlo 
(histogram). The samples are normalized to each other; this 
normalization is used for the final background estimates. 



the total number of data events with two or more tracks. 
Figure || shows a comparison of the data and MC dis- 
tributions for this sample. The preliminary error on 
this normalization is 9%. Two alternate normalizations 
were used as checks. The first normalized to protons 
on target using the decay channel mass distribution and 
(anti-) neutrino cross-sections. The second normalized 
to charged-current interactions in the calorimeter, and 
scaled by the decay-channel-to-calorimeter mass ratio. 
These cross-checks had uncertainties of 16% and 12% re- 
spectively, and were in agreement with the primary nor- 
malization. 

Monte Carlo events were compared to data as a check 
on the quality of the simulation. In such a comparison, 
the challenge is to isolate events of sufficiently similar 
topology to verify the Monte Carlo calculation of the 
background and at the same time maintain high statis- 
tics in the data sample. We used two methods to achieve 
high statistics comparison samples of events with two or 
more tracks. 

For the first sample, the vertex was required to be 
within the decay channel fiducial volume, with |a;| < 
127 cm, |y| < 127 cm, allowing the z position to be ei- 
ther in the chambers or the helium. Tight track angle 
cuts were imposed to remove cosmic rays, and a strict re- 
quirement on veto system activity was used to remove up- 
stream interactions. Other cuts on reconstruction, parti- 
cle identification, and vertex fit quality were removed. 
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Monte Carlo.) 
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FIG. 6. Top: longitudinal vertex position for events in the 
helium, where helium is defined as an event having a vertex 
more than 101.6 cm away from a chamber. (Crosses: data; 
histogram: Monte Carlo.) Bottom: ratio of data to Monte 
Carlo events. 

The majority of these events were from interactions in the 
chamber material or from interactions in the laboratory 
floor. In the data, from 502 events, 169 events had ver- 
tices reconstructed in the helium, defined as > 101.6 cm 
for the nearest drift chamber. This can be compared 
to Monte Carlo, which predicted (525 ± 84) events with 
(159 ± 25) events reconstructed in the helium. Com- 
parisons of data to Monte Carlo vertex distributions are 
shown in Figures 0, a, and pi 

The second sample was comprised of Monte Carlo 
events passing the five normalization cuts. Figure M 
shows a comparison of some event variables for this sam- 
ple. There is good, qualitative agreement between the 
data and MC in these and additional distributions. For 
events with a z vertex more than 101.6 cm from the cham- 
bers (helium events) the Monte Carlo predicts 28 events; 
40 data events are observed. 



TABLE IV. Estimated rates of background to the 
A'^" — > ^JL^iy) search 

Source l^fJ-i^) events 



DIS events 
DifFractive charm 
Difi'ractive n 
DifFractive K 
K1 decays from berm 
Other sources 



(3.9 ± 0.9) X 10-2 

(1.1 ± 0.1) X 10"^ 

(1.7 ± 0.1) X 10-* 

(3.3 ± 0.3) X 10-^ 

(3.9 ± 3.9) X 10-* 

< 2.5 X 10"* 



Total t-iniu) Background 



(4.0 ± 0.9) X 10-2 
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FIG. 7. Distributions comparing data (points) and Monte 
Carlo (histogram) for charged-current DIS interactions in the 
decay channel chambers. The distributions are: (a) transverse 
mass; (b) hadronic energy in the front of the calorimeter; (c) 
vertex pseudo-x'^/dof; (d) number of reconstructed tracks in 
the X view. 



After all cuts, the preliminary expected background is 
0.040 ± 0.009 events in jj/j, mode, 0.14 ± 0.02 events in ij,e 
mode and 0.13 ± 0.02 events in inr mode. As an example 
of the relative sizes of the contributions discussed above, 
the background sources for the N*^ —>■ ^ifJ.{v) mode are 
broken down in Table IV. 



VI. CROSS-CHECKS USING DATA 

Before looking at the data in the signal region, we per- 
formed a series of analyses on other fiducial and kine- 
matic ranges which gave us confidence in our Monte Carlo 
predictions. We point out that the two previous pub- 
lished analyses in the 0.3 to 3.0 GeV/c^ and low mass re- 
gion were examinations of other such kinematic regimes. 
Extra studies performed for this analysis included us- 
ing 1) identical analysis cuts applied to events within 
±15.2 cm of a chamber (the chamber region); 2) the 
chamber region with loosened cuts to increase inr ac- 
ceptance; 3) the "intermediate region" between 15.2 and 
101.6 cm from the chambers, with otherwise standard 
analysis cuts; and 4) events with well-reconstructed two- 
track vertices where the tracks were both identified as 
pions. The results were within 1.5cr of prediction in the 
above cases. 



TABLE V. Number of events in drift chambers which pass """■ siss Event: 3846 igate: i Date: Tue :ui 2212:27:121997 
Decay Mode 



Predicted Events 



Observed Events 



HH (chamber) 


1.6 


fie (chamber) 


1.8 


fj,TT (chamber) 


2.7 



We use the chamber region events as an example of 
these studies. This is a very powerful data sample be- 
cause if any observed events in the decay region were 
due to neutrino-He interactions, then there should be a 
factor of 28 more events in the chambers, after scaling for 
acceptance and mass. The numbers of observed events 
for the iifj,, /ie, and /itt modes are listed in Table ^ along 
with the prediction for neutrino deep-inelastic scatter- 
ing in the chambers. The observed events are consistent 
with the prediction, giving no indication of unexpected 
"clean" , two-track neutrino interactions in the chambers. 

VII. RESULTS OF THE SEARCH 

Using the signal event selection criteria given above, 
the Monte Carlo background predictions are given in Ta- 
ble VL The number of observed events is also shown. 



Three /^/i(j^) events were observed, which is considerably 
above the predicted background. No /le or /xtt events are 
observed, consistent with expectation. 

TABLE VI. Predicted and observed events passing all sig- 
nal cuts 



Decay Mode 



Predicted Events 



Observed Events 



titi{u) 


0.040 ± 0.009 


3 


fie{u) 


0.14 ±0.02 





flTT 


0.13 ±0.02 







FIG. 9. Run/Event 6133/3846: ///i(^) data event passing 
final cuts. 



Run : 6013 Event : 2 19S63 Igate 



Date: Fri May 2 09:' 




FIG. 10. Run/Event 6013/219863: fifi{iy) data event pass- 
ing final cuts. 



Run: 5S35 Event: S1705 Igate: 1 Date: IVed Dan 22 1S:23:07 1997 




FIG. 8. Run/Event 5835/81705: fifi{u) data event passing 
final cuts. 



The events are shown in Figs. WM A summary of 
the event rec onstr uction characteristics is shown in Ta- 
bles ^^ and [VIII|. Because of the large multi-hit TDC 



ring-buffer used for the veto wall, it is typical for events 
to have an average of 1.7 counters firing per neutrino in- 
teraction. For event 5835/81705, the TDC times with 
respect to the trigger are +404, +536 ns, where the posi- 
tive sign indicates the hits occurred after the trigger; for 
event 06133/03846, the hit occurred at +24 ns; and for 
event 6013/219863, hits occurred at -256, +320, +320, 
and +1192 ns. Given the 3.8 ns timing resolution of the 
veto counters, these hits are well out of time. 



TABLE VII. Kinematic and reconstruction quantities for the three candidate 
A*'" -^ MM('^) events. The sign on the muon energy refers to the charge of the muon 
(if measured) . 



Run/Event Mode 



Epl E;j2 -Prmiss minv rUT 

(GeV) (GeV) (GeV/c) (GeV/c^) (GeV/c^) 



5835/81705 


V 


-11.1 


±2.56 


2.42 


1.10 


5.08 


6133/3846 


V 


-92.0 


±5.85 


1.41 


0.88 


3.08 


6013/219863 


V 


±48.0 


±4.34 


2.07 


3.57 


4.66 



TABLE VIII. Vertex information for the observed \i[i events; V> — probability of 
an event having this pseudo-x^/dof or greater. 



Run/Event 



Vx 

(cm) 



(cm) 



|2^vcrtcx ^^chambc 

(cm) 



pseudo-Xvcrtcx/rfo/ 

iv>) 



5835/81705 


-46.7 


3.6 


196 


6.3/9 (62%) 


6133/3846 


46.5 


-38.6 


792 


166/17 (5%) 


6013/219863 


-59.2 


14.7 


186 


22.6/10 (21%) 



VIII. CROSS-CHECKS BY RELEASING CUTS 

The observation of the three ^^ data events prompted 
further tests comparing data to Monte Carlo predictions 
with reduced cuts. These studies provide cross-checks 
on whether the signal was manufactured by the cuts, 
whether the Monte Carlo accurately models the data just 
outside the cuts, and whether there is any indication of 
excess background with less restrictive requirements. 

The first study examined how data and Monte Carlo 
rates varied as the cuts were loosened gradually. The 
expected background and number of observed events at 
each step of the process is shown in Table IX . The first 



step was to remove the "clean cuts." At this point, one 
additional fj,fi event was observed. The total is written as 
3+1, explicitly separating off the three candidates. The 
topology of this additional event was in agreement with 
the topology for a DIS event (extra track segments, excess 
of neutral energy, and extra hits), and with the Monte 
Carlo estimate of 0.25 events. No /ie or /itt events were 
observed which is consistent with expectation. Thus, 
except for the excess three events, at this step there is 
agreement between data and Monte Carlo. Continuing 
to remove the cuts sequentially, the fiducial region was 
increased to include the chambers; finally, all energy and 
PID cuts were removed. In both of these steps, the data 
continued to be in agreement with the Monte Carlo, ex- 
cept for the three candidate events. 



A second study considered the effect of loosening the 
cuts on the expectation for a Monte Carlo which com- 
bined N^ decays in nnv mode with the standard back- 
ground Monte Carlo. This study normalized the N'^ de- 
cay Monte Carlo to 3 events for standard cuts and then 
observed the expectation with the cuts sequentially re- 
leased as shown in Table ^. There is agreement of the 
data with the combined model but the low statistics are 
not definitive. 

A third cross-check involved alternately releasing and 
then returning individual cuts and cut-pairs. For most 
cuts, no extra events entered the sample. When the re- 
quirement of "no extra track segments" was removed, one 
fiTT event entered the sample. When the "no extra clus- 
ters" cut was removed, one event identified as fj,e entered 
the sample. This is consistent with the Monte Carlo ex- 
pectation, where the extra cluster is due to photons from 
a 7r° and these cause the event to look electromagnetic. 
When both of these cuts are removed, two more events 
enter the sample. Their topology is consistent with DIS. 
In summary, this cross-check did not reveal events which 
are similar to the candidates. 



TABLE IX. Observed events and expected background as cuts are sequentially released. 



Sequential Change to Cuts 


Event Type 


Number of Data Events 


Monte Carlo Prediction 


All cuts 


MM 


3 


0.04 




/iTT 





0.14 




/ie 





0.13 


Remove the "Clean Cuts" 


fifi 


3 + 1 


0.25 




)in 





0.23 




fie 





0.70 


Include the chamber region 


fifi 


3 + 2 


1.3 




flTT 


7 


7.8 




fie 


5 


4.7 


Release energy & PID requirements 


lifi 


3 + 3 


2.3 




flTT 


10 


16.0 




fie 


10 


10.8 



TABLE X. Predicted and observed fifi events as cuts are sequentially released. Prediction shows the 
expected background and the combination of signal (from A''" decay Monte Carlo) and background. The 
signal is normalized to 3 events with all cuts. 



Event Type 


Predicted Background 


A'^-decay 


A^-decay + Background 


Data 


fifi with "standard cuts" 


0.04 


3.00 


3.04 


3 


releasing the "clean cuts" 


0.25 


5.30 


5.55 


4 


Including chamber regions 


1.3 


6.1 


7.5 


5 


removing E & PID requirements 


2.3 


6.1 


8.6 


6 



10 



IX, 



CONSIDERATIONS UNDER AN N° 
HYPOTHESIS 



5 GeV N° -^ finv 



In many ways, the three fiji events are consistent with 
a N^ decay hypothesis. The events pass the analysis cuts, 
where the background is estimated to be 0.04 events. As 
expected for a decay relative to an interaction hypoth- 
esis, all three events occur well within the fiducial vol- 
ume away from the chambers and are evenly distributed 
throughout the decay channel. The transverse mass, in- 
variant mass, and missing Pt are all consistent with a 5 
GeV/c^ N° decay (Fig. y). 

Unlike the background, in both the NHL and neu- 
tralino models, one would expect the /itt rate to be highly 
suppressed relative to leptonic decays. However, for a 5 
GeV/c^ NHL model, one would expect 1.4 times more 
jiev events [0. This is not inconsistent with the obser- 
vation of no /ie candidates, but neither docs it provide 
direct support. A neutralino model, on the other hand, 
can accommodate the observation of either only /i/i or a 
combination of /i/i and /le candidates by selecting appro- 
priate couplings. 

Globally, the events have one feature which is improb- 
able for an A''" decay hypothesis. Figure 12 shows the 
muon energy asymmetry {\Ei — E2\)/{Ei + E2) for the 
5 GeV/c^ N'^ simulation compared to the three data 
events. All three events have a muon energy asymme- 
try which is greater than 0.85. The probability that this 
occurs in a weak decay hypothesis pG| is less than 0.5% 
(including acceptance). 

Individually, there are also particular characteristics 
to note in two of the events. Event 06133/03846 has 
relatively poor vertex pseudo-x^/dof; the probability of 
an event having a pseudo-x^/dof greater than or equal to 
this is 5%. In event 6013/219863, three of the four hits in 
chamber DK4 immediately downstream of the vertex are 
missing. The probability of this occurring in a random 
event is less than 3 x 10^*^. DK4 had persistent high volt- 
age problems throughout the run and the absence of hits 
in this chamber may be due to this. However, we have 
examined the drift chamber readout string and found no 
indication that the chamber was misbehaving for this 
event and we have observed good tracks through this 
chamber for events proximate in time. Finally, we have 
discovered some evidence for a "coherent" inefficiency, a 
linked inefficiency between the x and y views. Such an 
inefficiency may occur at a level of (1.0±0.7)x 10~^ prob- 
ability, again making the missing of the three hits on the 
tracks highly unlikely. This event also has a third track 
attached to a downstream vertex. To interpret this event 
as an A^" decay, the added track must be interpreted as 
a delta ray. Such an interpretation is consistent with the 
calorimeter data. Event 5835/81705 does not suffer from 
reconstruction issues. 
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FIG. 11. Kinematic distributions (transverse mass, invari- 
ant mass and missing transverse momentum) for the 5.0 
GeV/c^ N° Monte Carlo. The histograms show the MC; the 
arrows indicate the three observed events. 
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: Asymmetry = 
abs(E1-E2)/(EH-E2) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Muon Energy Asymmetry 

FIG. 12. The muon energy asymmetry 

{\Ei- E2\)/{Ei+E2). The histogram shows the 5 GeV/c^ 
a" MC;the arrows indicate the three observed events. 
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NuTeV limit at 90% C L. (Preliminary) 



00 




NuTeV Limit (B and Berm) (solid) 
DELPHI Limit (dashed) 



Mass of NHL (GeV) 



FIG. 13. NuTeV limit on NHL production Ironi B decays 
and berm production. 




decay length (km) 

FIG. 14. NuTeV limit on neutralino production. This limit 
is generic for an A'^" produced at the target. The right axis 
(# of events) follows the shaded contours. 

Interesting limits on the production of NHLs and neu- 
tralinos can be set that indicate the sensitivity of the 
decay channel search; these limits are set by calculat- 
ing one-sided limits using a frequentist approach with- 
out background subtraction. Because the source of the 
events is unclear, the unified approach of Feldman and 
Cousins [pl| has deliberately not been used. The limit 
on NHLs reaches mixing parameter values below |C/p = 
10~^ and is consistent with the results from the Delphi 
experiment at LEP |^ , as shown on Figure |3[ NuTeV 
is the first experiment to set limits on the production 
of long-lived neutralinos in this mass range which deca' 
by i?-parity violation. The result is shown in Figure 



cay 



This limit, although motivated by a neutralino hypothe- 
sis, is a generic limit applicable for any model of neutral 
particle production at the target B. 



X. CONSIDERATIONS UNDER VARIOUS NULL 
HYPOTHESES 

Several aspects of the candidate events are similar to 
those from neutrino interaction backgrounds, and might 
be indicative of unaccounted-for sources. First, all three 
events occurred during the high rate :/-mode as op- 
posed to F-mode running periods; the i^ to F event ra- 
tio is expected to be 4:1 for all events and 1.5:1 for 
low-multiplicity events. Second, if the events were pro- 
duced by neutrino interactions, then one would expect 
a high energy leading fi~ for ly-mode running. In the 
two cases where the charge of the leading muon can be 
measured, the charge is negative. The event kinemat- 
ics {Mt, Mf^f^, Pt) are consistent with the background 
simulation (Fig. |l5|) and, for DIS background, the prob- 
ability for three events with the observed energy asym- 
metry is 25 to 35% (Fig. |l^). If these are interpreted 
as DIS events, then the visible (measured) squared four- 
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FIG. 15. Kinematic distributions (transverse mass, invari- 
ant mass and missing transverse momentum) for fifi events 
from the background Monte Carlo. The histograms shows 
the MC while the arrows indicate the three observed events. 
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Muon Energy Asymmetry 

FIG. 16. The muon energy asym- 

metry {\Ei — E2\)/{Ei + £2). The histogram shows the fifi 
background Monte Carlo while the arrows indicate the three 
observed events. 

On the other hand, it is difficult to explain these events 
as neutrino interactions. If all of these events are due 
to charged-current neutrino interactions in the helium 
which give a prompt second muon (charm production, 
trident production, etc.), then the rate is two orders of 
magnitude greater than predicted by the Monte Carlo. 
If such a source was not included in our simulations, 
then we would also expect between 10 and 300 additional 
events at 99% confidence level pl| in the chambers and 
between 30,000 and 220,000 events at 99% confidence 
level in the calorimeter itself. We actually see only 2 fifi 
chamber events in agreement with our Monte Carlo (see 
Table [X) and 54 events with similar kinematics were 
found in the calorimeter [[l9| . In order to show the dif- 
ficulty of developing reasonable "null" (background) hy- 
potheses, we consider three examples below. 

The first null hypothesis is that the events were ac- 
tually prompt fj,fj, events produced in the chambers and 
misreconstructed in the helium. If we try to force ver- 
tex reconstruction for the two muons in 6013/219863 at 
the nearest downstream chamber (DK4), the pseudo- 
Xvertex/dof gocs from 22.9/10 for the standard fit up 
to 290.2/11 for the forced fit. If the third track is also 
forced to be on the vertex, the pseudo-x^gj.tgx/dof rises 
to 323.6/17. If one does not require the track to be in the 
chambers, but requires all three tracks to come from the 
same vertex, then one obtains 217.1/16. Likewise, event 
6133/3846 produces a poor fit to the chambers. For the 
chamber downstream of the vertex (DK5), the pseudo- 
Xvcrtcx/dof is 94,720/18 and for the testbeam chamber 
upstream it is 543/18. These are well above the origi- 
nal fit of 166/17. Finally, for 5835/81705, the standard 
fit gives 6.3/9 which increases to 51.5/10 if forced to be 
reconstructed in the chamber upstream (DK4). Aside 
from the problems of forcing a fit to the chambers, this 
hypothesis leaves unanswered the questions: 1) why are 
only /i^ events misreconstructed, and not fin or fxel and 
2) what is the source of these excess /z/i's? 

A second null hypothesis rejects the idea that the iifi 



events are prompt, and instead attributes them to unsim- 
ulated ly + X^fi + n + Y interactions, where the pion 
decays to produce the second muon and the other parti- 
cles in the reaction (signified by Y) are not seen. This 
has the advantage of being more plausible than an un- 
expected source of prompt fifi events, but it suffers from 
several problems. First, if these are events which occur in 
the chambers, then one must explain why the fi/j, events 
are misreconstructed in the helium. The transverse mo- 
mentum associated with pion decay is very small, so the 
muon track would generally point back to the original 
chamber vertex. For those events which do show an off- 
set, there will be a much higher probability of occurring 
in the 15.2 to 101.6 cm range than well beyond 101.6 cm, 
where these events are found. Also, only 7% of the pions 
will decay within the decay channel so many undecayed 
/ITT events should be observed. 

Quantitatively, the data in these other channels can be 
compared to different scalings of the observed three fifj, 
events. The DIS Monte Carlo predicts a ratio of fiir to fifx 
events of 3.25 for the sample with all the standard cuts. 
With this ratio, the probability of seeing 3 fifi events and, 
as observed, /xtt is 0.31%. Furthermore, one can scale 
these fifi events in the helium by the ratio of masses and 
acceptance for the chambers versus helium and compare 
to the observed /xtt and fifi events seen in chambers with 
the loosest cuts. Given that ratio, the 10 fin (2 fifi) events 
observed in the chambers are consistent at a confidence 
level of only 2.1 x lO"*^ (2.12 x lO'^) with the hypoth- 
esis that the three fifi events in the helium are due to 
an unaccounted-for fin source. It is therefore extremely 
unlikely that an unsimulated fin source is responsible for 
the three fifi events in the helium. 

A third null hypothesis is that the candidates are neu- 
trino interactions in the chambers producing a fiK fi- 
nal state followed hy a K ^ fiiy decay that causes the 
event vertex to be misreconstructed downstream in the 
helium region. This process has the advantage that the 
decay angle is sufficiently large for an event produced in 
the chambers to reconstruct in the helium. Fits to the 
three candidate events under this hypothesis show that 
that the fits are kinematically possible. Taking the kaon 
momenta calculated from the fit, the probabilities are 
17.6%, 11.6% and 8.2% for the kaon to decay. Because 
the decay channel tracking is largely located downstream 
of the decay region, most two-track vertices, including 
most simulated N^ events, can also be reconstructed as 
displaced-vertex kaon decays with high quality. 

However, the kaon hypothesis does not explain these 
events for several reasons. First, direct production is 
Cabibbo-suppressed, so most production is through frag- 
mentation. Direct production is also most likely in D 
running, whereas the events are only seen in v running 
periods. Approximately 20% of events produce kaons 
through fragmentation, but this occurs mainly in high 
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multiplicity events which would be cut by the "clean 
cuts." This process is included in the 0.040 event back- 
ground estimate, to which it contributes ~0.0025 events. 
Clean events could be produced by difFractive produc- 
tion; the large angle of the high energy muons in these 
events, however, is uncharacteristic of diffractive produc- 
tion and the excess is far above expectation (see Table 

0)-. . . 

Finally, as with the pion decay hypothesis only a small 

fraction, 22%, of the kaons from neutrino interactions 
will decay. In addition, because kaon decays have large 
Pt, only 45% of the fi — K producing two muons in 
the calorimeter will have a vertex formed by the two 
muons. In the other cases, the decay plane is aligned 
such that the two muons fail to verticize in the decay 
channel. (Such unaligned dimuon events consistent with 
all but the vertex cut are not observed in the data.) As 
in the case of the pions, we can use scaling arguments 
to compare the observed data to the predictions of this 
hypothesis. In the data, there are one /Lte and two fin 
chamber events with all the standard cuts but with the 
track energy requirements removed. The scaling factor 
is 8.24 between fiK events in the drift chambers and fifi 
events in the helium. Under this hypothesis, the proba- 
bility to see, as observed, three or fewer /x-hadron events 
in the chambers given the three /z^ events in the helium is 
0.37%. Therefore, the paucity of /z-hadron events in the 
chambers excludes at 99.6% confidence level the possibil- 
ity of an unsimulated fiK source large enough to explain 
the three observed fifj, events in the helium. 



XI. SUMMARY AND CONCLUSIONS 

In summary, NuTeV has observed 3 fj-ii events, ijtt 
and fie events with transverse mass above 2.2 GeV/c^. 
The expected backgrounds were 0.040±0.009, 0.14±0.02, 
and 0.13 ± 0.02 events (preliminary), respectively. 

NuTeV has set new limits on NHL and neutralino de- 
cays based on this analysis. The NHL limits are consis- 
tent with the Delphi result. The neutralino results are 
the first in this kinematic region for the long-lived x'^ 
which decays with i?-parity violation. 

In conclusion, the rate corresponding to the observed 
three events is not consistent with Standard Model pro- 
cesses we have identified and the source of the events is 
not clear. 
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